Time and length scales in supercooled liquids 
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We numerically obtain the first quantitative demonstration that development of spatial correla- 
tions of mobility as temperature is lowered is responsible for the "decoupling" of transport properties 
of supercooled liquids. This result further demonstrates the necessity of a spatial description of the 
glass formation and therefore seriously challenges a number of popular alternative theoretical de- 
scriptions. 
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Transport coefficients in liquids approaching the calori- 
metric glass transition change by many orders of magni- 
tude [l| . It was discovered about a decade ago that con- 
ventional hydrodynamic relations are not valid in super- 
cooled liquids in the sense that various transport prop- 
erties "decouple" from one another 0,0. In a stan- 
dard liquid the viscosity, 77, translational diffusion coef- 
ficient, D, and temperature, T, are linked by the hy- 
drodynamic Stokes-Einstcin relation, Drj oc T, which 
breaks down for supercooled liquids. The product Di] 
can be orders of magnitude larger than its hydrodynamic 
expectation at the glass transition, a puzzling observa- 
tion which has received considerable interest in the last 
decade 0, H 0, H, H . 

We have numerically obtained the first quantitative 
link between this decoupling phenomenon and the ex- 
istence of spatial correlations in the dynamics of super- 
cooled liquids, also called "dynamic heterogeneity" Q. 
We define here dynamic heterogeneity as expressing the 
fact that local dynamics of the liquid becomes spatially 
more correlated as T decreases |7|. This naturally im- 
plies the existence of a wide distribution of relaxational 
time and length scales 0, Q . Decoupling results because 
different transport coefficients correspond to different av- 
erages over those broad distributions [5|. 

Increasing heterogeneity was previously reported in a 
number of numerical works but its direct relevance 
to transport was not established. On the other hand, ex- 
periments have quantified decoupling 0, 0, 0] vi 311 * s P a " 
tial correlations were only indirectly measured 0, [H| be- 
cause they are presently inaccessible to scattering ex- 
periments since neutrons/light probe too small/large 
length scales, and no obvious temperature dependence 
was found 0, |n| . The novelty of this work is therefore 
to reconcile in a direct manner two facets of the slow dy- 
namics of supercooled liquids. Our results demonstrate 
the correctness of the widely shared belief that decou- 
pling is most naturally interpreted in terms of heteroge- 
neous dynamics, a result which has been actively seeked 
in the last decade and has deep theoretical implications. 

Our argument shall be presented in three steps. In 
large scale numerical simulations of a well-characterized 
liquid model [l2T ]. we first quantify dynamic heterogene- 



ity, then decoupling, to finally establish in a third step the 
quantitative connection between these two aspects. We 
finally discuss the important theoretical consequences of 
our findings. 

We investi gate the binary Lennard- Jones mixture pro- 
posed in Ref. |l2j with Na — 1097 particles of type A and 
and Nb — 275 particles of type B at density p = 1.2. The 
N = Na + Nb particles interact via a Lennard- Jones po- 
tential V{v a p) = 4e Q/3 [(<J a p/r a p) 12 - (a- al3 /r a p) 6 ] , with 
a, (3 = A,B. Time, energy and length are measured 
in units of cfaa and caa, and \JmA<J 2 A A l £aa, respec- 
tively. Other parameters are €ab = 1-5, £bb = 0.5, 
&bb = 0.88, a ab — 0.8. Newton equations are inte- 
grated via a leapfrog algorithm with time step 0.01. Ve- 
locity rescaling is used to thermalize the system. We 
study the system at equilibrium for a wide range of tem- 
peratures, T 6 [0.42,2.0]. Since we measure fluctuations 
of local dynamical quantities, extremely long simulations 
are necessary to ensure not only thermal equilibrium, but 
also that sufficient statistics is recorded. We use the par- 
allelized algorithm developed by Plimpton 0], and at 
all T runs of length at least 100r Q are performed. Our 
computer capabilities fix the lowest studied temperature. 
Characteristic temperatures for this system are the on- 
set of slow dynamics, T a w 1.0, and T c w 0.435, the 
location of the mode-coupling singularity in the analysis 
of Ref. [13 . 

First, we measure the temperature dependent coher- 
ence length, £(T), associated with the ordering of the 
liquid's dynamics by measuring spatial correlations be- 
tween individual particle relaxations [9j. We use 
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6^(0) - v) cos (tk-[r J (t)-v J (0)}) (1) 



as a natural local indicator of the dynamics, since 
F s (k, t) = (Fk(r, t)) is the real part of the standard self- 
intermediate scattering function; (■ • • } stands for an en- 
semble average at temperature T, while rj(t) is the po- 
sition of particle j at time t. In our definition, dynamic 
heterogeneity implies that Fy(r, t) becomes long-ranged 
correlated as T is lowered [14j , as can clearly be observed 
in the snapshots of Fig. ^ 
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FIG. 2: Dependence of the coherence length £(T) measuring 
the mean size of the clusters of Fig. Qon the inverse temper- 
ature from the spatial fluctuations of F^(r, t), for k = 5.41 
and t = t, see Eq. @. These data are discussed in detail in 
Ref. El. 





FIG. 1: Snapshots of the simulation box where points with 
SF — i<k(r, t) — F s (U.,t) > are represented with radii pro- 
portional to SF at temperatures T = 2.0, 0.6 and 0.45 (top to 
bottom), with k — 5.41, t = r(fc = 5.41, T). Increasing spa- 
tial correlations of the particles' individual dynamics when T 
is lowered is evident. 



The measurement of the mean size of the dynamic clus- 
ters shown in Fig. ^ involves the study of a two-point, 
two-time correlation function as already discussed in sev- 
eral papers 0H El El El We extract 1{T) from the 
wavevector dependence of the Fourier transform of the 
correlator Ck(r): 



Ck(r) 



(F k (0,r)F k (r,r))-(F k (r,r))- 
(F k (r,r)2)-(^ k (r,r))2 



(2) 



where r = r(fc = k|,T) is the relaxation time defined 
in a standard way from the time decay of F s (k, t) E3 • 
The shape of C k (r), the temperature dependence of t(T), 
and their theoretical interpretation in the context of a 
renormalization group analysis are the main object of 
Ref. El, 

so that we only report the data for t(T) in Fig. [21 
without discussing them further in the present report. 
Note that this first step is also the most demanding in 
terms of numerical resources. 

Secondly, to probe decoupling, we measure with a great 
precision the temperature and wavevector dependences 
of the relaxation time r(k, T) defined above. Our results 
are presented in Fig. We find that the temperature 
dependence of r(fc — * 0, T) is the same as the inverse 
diffusion constant D -1 , as expected in a diffusive regime, 
while for wavevectors close to the first peak of the static 
structure factor, kg rj 7.2, the temperature dependence 
is stronger, and follows that of the viscosity, establishing 
decoupling in our model system. This finding has alr eady 
been reported in several numerical works 

U II El El, 

and is not unexpected in a system characterized by broad 
distributions of time scales. 

We can now establish our main result which is the link 
between decoupling and dynamic heterogeneity. It sterns 
from intermediate wavevectors, < k < ko, for which 
a crossover is observed in the temperature evolution of 
r(fc,T), uniquely governed by the value of the scaling 
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FIG. 3: Temperature dependence of r(k,T) for various 
wavevectors, k = 7.21, 6.61, 6.00, 5.41, 4.81, 4.21, 3.00, 2.40, 
1.80, 1.20, and 0.60 (bottom to top). The dashed line is 10/D, 
indicative of the r(k — > 0, T) limit. Decoupling is observed as 
the temperature dependence is stronger at larger wavevectors. 
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FIG. 4: The quantity X(k,T), Eq. ©, as a function of the 
scaling variable k£(T) for various T. The horizontal full line 
is the diffusive prediction, rDk 2 ~ const. Departure from 
diffusion arises at large kl, directly demonstrating that the 
decoupling of different transport properties results from an 
increasingly spatially correlated dynamics. 



variable k£(T). We find that r(fc,T) follows D^ 1 when 
k£ < 1, or 77 for k£ > 1. It is therefore useful to define 
the following ratio, 



X(k,T) 



r(k,T)D(T) 
r{k,T )D(T o y 



(3) 



since X(k, T) is wavevector independent through the de- 
nominator, and is temperature independent if r(fc, T) oc 
D~ l (T), so that by definition X(k,T) = 1 if no decou- 
pling occurs. 

Instead we find that time scales r(k, T) spanning 7 
orders of magnitude can be collapsed on a unique, non- 
trivial curve, as shown in Fig. so that 



X(k,T) ~ X[k£(T)]. 



(4) 



We find that X(x) = 1 + x 13 with f3 w 1.6 represents the 
numerical data quite well, so that the diffusive regime is 
limited to small values of the scaling variable kl. Note 
that this scaling extends from the onset temperature T a 
to well below the mode-coupling temperature T c , From 
the scaling behaviour Q, our main conclusion is there- 
fore that dynamic heterogeneity emerging at T a and in- 
creasing when T is lowered is directly responsible for de- 
coupling. 

Decoupling phenomena have been experimentally char- 
acterized close to the calorimetric glass transition where 
they are more pronounced, as can be understood from 
Fig. We believe, however, that the behaviour is qual- 
itatively similar to our numerical investigations because 
even an increase of several decades in time scales corre- 
sponds to a very modest change of the coherence length 
£(T), which is the key quantity of the problem. Exper- 
iments report moreover that Drj ~ r/ a with 1 > a > 
00,13 We find similarly that, at fixed k, tD ~ $ . 



The dynamic scaling discussed in Ref. [lif. 
yields indeed the observed power law, 



Dr/ ~ i] 



-l/z 



(5) 



with an exponent f3/z rs 0.35. Various values for this ex- 
ponent have been reported from experiments [illl8| and a 
precise characterization of all these exponents for various 
liquids on a wide range of length scales would be most 
useful. Note that the use of the crossover scaling func- 
tion X{x) in Eq. JSJ) would automatically yield a smaller 
effective value of the exponent 0/ z on a restricted time 
window. 

Finally, we discuss the theoretical interpretations and 
consequences of Fig. 0] These results constitute the first 
direct demonstration that, much as in ordinary critical 
phenomena, large spatial correlations not only accom- 
pany but also influence the properties of supercooled liq- 
uids, and therefore the glass formation itself. Thus, they 
are a striking confirmation that spatial approaches are 
necessary to understand the formation of glasses, and 
several such quantitative app roaches can be found in the 
literature [H E |U |U H|| . 

When a growing length scale for dynamic heterohene- 
ity is included, our results are very naturally explained, 
just as the decoupling phenomenon is 0. Tuning the 
wavector in Q amounts to probing the dynamics on dif- 
ferent length scales, smoothly interpolating between Z> _1 
(k — * 0) and r] (k » k ), as is observed in Fig. 0] which 
clearly demonstrates that the crossover, k£(T) ~ 1, is 
ruled by the increasing correlation length of dynamic het- 
erogeneity. 

Our findings also confirm that dynamic heterogeneity 
is a central aspect of the dynamics of supercooled liq- 
uids in that time and length scales are intimately con- 
nected 0, El ■ This is at odds with the opposite belief 
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that slow dynamics emerges because of the local blocking 
of the particles without any relevant length scale beyond 
the interparticle distance, the famous "cage effect" [23j. 
Our results make it clear that the alpha-relaxation is in- 
stead a cooperative phenomenon where single particle dy- 
namics are coherent on the length scale £(T) much larger 
than , which is in turn directly responsible for the 
temperature behaviour of the time scales. Moreover, l(T) 
starts to grow significantly and connects to time scales 
even in the regime T c < T < T a where a mode-coupling 
analysis supposedly applies 0, Hj| . This confirms that 
heterogeneous dynamics, decoupling and activated dy- 
namics set in at the onset temperature T D which is there- 
fore the key temperature scale of the problem, as opposed 
to T c where no significant change of mechanism takes 
place 0, 01 . We note that the absence of non-trivial 
spatial correlations [25| , the incorrect identification of T c 
as a key temperature scale, and the absence of decoupling 
in the present formulation of the mode-coupling theory 
represent major failures of this approach. 

Our results constitute therefore a sharp new test to 
discriminate between the many theoretical approaches to 
the glass transition problem lj. Indeed, any theory in 
which time scales do not directly follow from the exis- 
tence of spatial correlations growing when T is decreased 
below the onset temperature T Q is seriously challenged 
by the present work. Quite importantly, this includes a 
number of approaches which have been nonetheless much 
applied to describe experimental results |U |2f| |2]j ■ 

Our conclusions are drawn from a scaling relation dis- 
covered for a specific, yet paradigmatic, model system. 
Although we believe they are generic to supercooled liq- 
uids, the absence of any data on different systems em- 
phasizes both the novelty of this work and the need for 
further detailed investigations of the slow dynamics of 
various supercooled liquids. It is clear, for instance, that 
the precise experimental characterization of £(T) should 
become a central goal for future investigations, and the 
present work therefore suggests a new way to access phys- 
ically relevant length scales. 
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